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Abstract 

The  potential  wave  propagation  is  described  by  defining  the  driving  parameters  of  the  propagation 
characteristics.  Then  a parametric  study  with  a non  linear  code  is  performed  in  order  to  define  the  mesh 
refinement  needed  to  obtain  a sufficient  accuracy  in  the  predictions  of  inflow  and  outflow  perturbation 
propagation.  The  results  are  compared  with  a linear  version  of  the  code  as  well  as  with  a linear  code  for 
flat  plate  response  showing  the  good  behaviour  of  the  non  linear  code.  As  additional  test  cases,  the 
behaviour  of  the  code  for  simulating  the  response  o a cascade  of  flat  plates  to  potential  disturbances  is 
investigated.  Results  from  the  different  linear  and  non  linear  codes  are  also  compared.  Different  cases  are 
investigated  including  one  with  large  resonances,  which  is  shown  to  be  the  most  difficult  case  to  be 
predicted,  in  special  the  resonance  peaks  as  well  as  the  trailing  edge  region  regarding  the  wave  phase 
behaviour. 


1 Introduction 

The  principal  aircraft  and  engine  manufacturers  in  Europe  are  facing  increasing  pressure  to  reduce  aircraft 
noise  levels.  Regarding  the  noise  from  engines,  low  pressure  turbine  (LPT)  noise  is  becoming  more 
important  due  to  the  decrease  of  jet  noise  with  the  high  by  pass  ratio  cycles  and  lower  noise  from  the 
advanced  design  fans  and  negative  scarfed  intake  angles.  In  a LPT,  the  interaction  of  unsteady  flow  with 
the  turbine  blading  produces  fluctuating  forces  along  the  blades  and  results  in  unwanted  effects  such  as 
noise  and  forced  vibration.  Such  excitations  are  mainly  generated  at  multiples  of  the  blade  passing 
frequency  and  arise  from  a variety  of  sources.  The  two  principal  types  of  such  interaction  are  usually 
referred  to  as  potential  flow  and  wake  interaction.  The  former  is  associated  with  unsteady  pressure 
variations  generated  by  the  upstream  and  downstream  rows  and  could  be  of  serious  concern  when  the 
axial  spacing  between  neighbouring  blade  rows  are  small  or  flow  Mach  numbers  or  excitation  frequencies 
are  high.  The  wake  interaction  is  due  to  vortical  perturbations  generated  by  wakes  convected  from  the 
viscous  flow  on  upstream  rows. 

Trying  to  understand  the  generation  of  noise  by  the  mechanisms  above  presented,  simulation  of  unsteady 
aerodynamics  with  linear  and  nonlinear  CFD  codes  is  an  ongoing  activity  within  the  turbomachinery 
industry.  However,  the  accurate  resolution  of  the  unsteady  aerodynamics  and  acoustics  on  low  pressure 
turbines  remains  a formidable  computational  problem.  This  is  due  mainly  to  the  spatial  resolution 
required  to  accurately  simulate  the  incoming  waves  (potential  and  vortical  fluctuations)  and  the  radiated 
cut-off  and  cut-on  waves,  even  in  the  case  that  only  the  2D  problem  is  considered.  Traditionally  linear 
codes  have  been  used  in  the  simulation  of  unsteady  flows  in  turbomachinery.  This  approach  allows  the 
use  of  relatively  coarse  meshes  and  they  need  less  computational  effort  when  compared  with  non  linear 
Euler  solvers.  However  they  assume  that  perturbations  are  small  when  compared  with  steady  variables 
and  they  need  to  solve  a new  set  of  equations  for  the  unsteady  variables.  This  approach  has  been  taken  by 
many  authors  ([  1 ] and  [ 2 ] among  others). 

A different  approach  is  taken  at  ITP,  where  a non  linear  unstructured  Navier  Stokes  solver  is  intended  to 
be  used  for  resolution  of  both  the  steady  and  unsteady  aerodynamics  and  acoustics,  that  will  allow  the  use 
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of  only  one  code  to  solve  both  problems  coupled.  This  seems  to  be  the  current  approach  taken  by  some 
other  authors  to  simulate  wake-boundary  layer  interaction  [ 3 ] and  to  simulate  potential  interaction  [ 4 ]. 

For  this  investigation  an  in-house  code,  with  the  non  linear  MU2S2T  and  linear  MU2S2TL  versions  is 
used.  An  extensive  investigation  on  the  different  parameters  involved  in  the  propagation  of  waves  is 
being  carried  out  with  the  objective  to  assess  the  mesh  resolution  requirement  and  the  capability  of  the 
solver  to  accurately  simulate  wave  propagation  and  unsteady  response  of  blade  rows. 

The  results  will  be  presented  in  an  step  by  step  approach,  from  the  different  parameters  involved  in  noise 
propagation  to  the  comparison  with  analytical  linear  solution  of  a flat  plate  response  to  incoming 
disturbances. 


2 Theoretical  background  of  the  wave  propagation 


The  equation  governing  the  unsteady  velocity  potential  for  small  deviations  for  an  uniform,  isentropic  and 
irrotational  mean  flow  has  the  following  expression. 


d20  . M 320  „ 2 x 

— _ — 2 + (1-A/  ) — r- 

dr 2 c dxdy  dx2 


= 0(1) 


that  has  expression  of  the  type 


<t>  = YZAnje(COn,*kyjy+kxn'iX)  ( 2 ) 

as  solution. 

For  turbomachinery  kind  of  flow  fields,  the  temporal  and  circumferential  eigenvalues  are 

o)n  = nCl  ( 3 ) 


Ay 


and  a expression  for  the  axial  wave  number 


0, 


Af  x (kyjM  y +-)± 


(kyJMy  + ®)-ky(l-Ml) 


1-M2 


(5) 


From  (5),  the  resonance  frequency  can  be  defined  as 


' ckyj(l-Mz) 

°>r=- 7-2 — (6) 


defining  the  cut-off  ratio  as 


f = — (7) 
cor 


Ratios  higher  that  unity  give,  according  to  (5),  the  so  call  cut-on  waves  (waves  propagating  with  constant 
amplitude),  while  ratios  lower  than  unity  give  cut-off  waves  (exponentially  decaying  waves  in  axial 
direction). 


For  cut-on  waves,  the  wavelength  in  axial  direction  can  be  defined 
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for  a propagation  wave  number 


3 Wave  propagation  simulation  with  CFD 

From  the  expressions  derived  above,  it  can  be  easily  seen  that  the  physical  characteristics  of  the 
propagation  of  a potential  wave  depends  on  the  following  parameters. 

O = / (Ux , Uy , p, T, Rg ,y, M, a, Ay)(  10  ) 

Wliere  y stands  for  any  unsteady  fluid  variable. 

That  involves  longitudes,  time,  mass,  temperature  and  angular  quantities.  If  the  speed  of  sound  c,  Rg,  the 
wave  length  Xy  and  the  density  p are  taken  for  nondimensionalization,  the  problem  is  reduced  to 

= Re)  (11) 

Mx  and  My  are  mach  number  in  x and  y direction,  % is  the  cut-off  ratio,  or  the  ratio  of  axial  to 
circumferential  wavelength  (that  involves  also  the  Mach  numbers,  but  gives  the  axial  propagation 
characteristics),  and  Re  states  for  the  Reynolds  number  based  in  the  wavelength  Xy. 

If  the  propagation  is  to  be  simulated  with  a CFD  code,  parameters  implicit  to  the  code  will  be  also 
relevant  to  the  problem.  These  parameters  have  to  do  with  the  time  and  space  discretization  as  well  as  the 
scheme  used  for  the  stabilization  of  the  system  of  equations. 

0 = f(Mx,My,^,Re,CFL,a2,cr4)  (12) 

Other  set  of  characteristics  in  the  simulation  of  the  wave  propagation  are  the  parameters  depending  on  the 
definition  of  the  mesh  employed  for  the  computation,  typically  Ny  cells  in  y direction  and  Nx  cells  in 
axial  direction.  For  wave  propagation,  is  of  interest  to  define  ny  and  nx  as  the  points  per  circumferential 
and  axial  wavelength  respectively. 

<t>  = f(Mx,My,^,Re,CFL,cr2,a4,nx,nv)  (13) 

In  real  turbomachinery  problems,  potential  and  vortical  wakes  perturbations  are  present.  As  this  paper  is 
focused  on  potential  disturbances,  that  are  not  affected  by  definition  by  the  Reynolds  number,  no  analysis 
on  this  parameter  will  be  done  here.  Regarding  the  effect  of  the  artificial  viscosity  parameters  included  in 
the  code  for  stabilization,  namely  02  and  a4,  their  effect  will  not  be  investigated  since  has  been  already 
investigated  on  the  same  code  by  Corral  et  al.  [ 6 ] . Nevertheless  some  conclusions  of  their  work  will  be 
used  and  commented  here.  The  only  effect  of  the  CFL  number  on  wave  propagation  is  computing  time 
and  no  further  results  on  the  simulation  will  be  presented  here. 

Therefore  the  parametric  investigation  will  cover  the  dependence  of  the  results  on  the  following 
parameters: 

® = f(Mx,M y,^nx,ny)  ( 14) 

In  addition,  studies  on  the  effect  of  the  quasi-structured  versus  fully  unstructured  meshes  has  been 
performed.  For  this  analysis,  quasi-structured  meshes  composed  of  quadrilateral  cells  with  one  preferred 
diagonal  direction  have  been  used. 


In  the  Figure  1,  the  results  for  a potential  perturbation  imposed  in  the  outlet,  with  a ratio  of  5.0,  M=0.4 
are  displayed.  Due  to  the  high  propagation  angle  of  the  wave,  the  simulation  presents  a low  number  of 
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points  per  axial  wavelength  and  the  results  present  a deficient  prediction  of  the  cut-on  wave  (the  mesh  has 
80  points  in  circumferential  direction  and  approximately  10  points  per  axial  wavelength).  Furthermore  the 
Figure  1 shows  that  there  is  some  effect  of  the  chosen  diagonal  direction  on  the  damping  of  the 
propagating  wave. 

ffiw-e  hropngaHon  Effect  of  ihf  meah. 


Figure  1:  Effect  of  the  mesh,  f:  5.0,  0):  195793.4  rad/sec,  outlet  perturbation 

In  turbomachinery  noise,  several  circumferential  modes  are  usually  propagated,  and  therefore  the  effect  of 
a preferred  diagonal  direction  would  have  a different  effect  on  each  mode,  hence  worsening  the  accuracy 
of  the  simulation  for  some  of  the  modes.  In  order  to  avoid  that,  the  use  of  fully  unstructured  meshes  is 
recommended,  trying  to  minimise  the  numerical  dissipation  of  the  propagation  waves  due  to  cell 
structure. 

All  the  computations  shown  in  this  study  have  been  performed  using  the  scalar  dissipation  model. 
Therefore  limiting  the  test  cases  to  a range  on  Mach  number  of  0.3<M<0.7.  However  when  having  very 
low  or  high  subsonic  Mach  numbers  the  matricial  dissipation  model  should  be  used,  that  has  been  already 
shown  by  to  significantly  improve  the  simulation  of  propagation  of  entropy  and  pressure  waves  which 
are  more  damped  when  using  a more  classical  scalar  model. 

One  of  the  critical  features  involved  in  the  simulation  of  wave  propagation  with  CFD  is  the  numerical 
boundary  conditions  implemented.  As  recognised  by  other  authors  [ 1 ],  the  ID  unsteady  non  reflecting 
boundary  conditions  (UNRBC)  are  insufficient  for  an  accurate  simulation  of  propagating  waves  in  short 
domains,  and  therefore  more  elaborated  boundary  conditions  are  needed.  In  this  investigation,  the  ID  and 
2D  UNRBC  by  Giles  [ 7 ] have  been  implemented  and  their  influence  on  the  wave  propagation  is  shown. 


Effect  of  the  Unsteady  Non  Refecting  Boundary  Conditions 


Figure  2:  Effect  of  UNRBC  with  £ 1.1,  M:  0.4 
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Figure  2 shows  that  in  the  case  of  the  ID  non  reflecting  boundary  conditions,  the  reflection  can  be  up  to 
60  % of  the  downstream  wave,  giving  a constructive  interference  and  a pressure  field  with  peak  to  mean 
values  close  to  50  % higher.  In  all  the  analysis  presented  hereafter,  2D  non  reflecting  boundary  conditions 
will  be  used. 

The  effect  of  the  cut-off  ratio  B,  has  been  investigated  through  the  simulation  of  different  cut-off  and  cut- 
on  waves.  To  avoid  crossed  effects  with  the  mesh  definition,  all  the  simulations  have  been  carried  out 
with  the  same  mesh. 

For  the  postprocessing,  the  wave  splitting  technique  proposed  by  Wilson  [ 8 ],  has  been  implemented, 
which  allows  to  investigate  in  detail  the  behaviour  of  downstream  and  upstream  waves,  by  means  of 
eigenvalues  and  eigenfunctions. 

Figure  3 shows  results  of  some  simulations  of  waves  with  E,=0.7,  1.1,  3.0  and  5.0  imposed  both  in  the 
inlet  an  in  the  outlet,  with  ny=80  points  and  Nx  constant  to  320  points 


WAVE  PROPAGATION.  CirzvnferenUal  points  per  uiovtlenghl  80. 


Figure  3:  Wave  propagation  for  £,:  0.7,  1.1,  3.0,  5.0  and  ny:  80  points.  Nx:  320 points 

As  it  can  be  seen,  the  waves  with  frequency  ratios  of  0.7,  1.1  and  3.0  are  in  general  well  simulated,  with 
an  small  axial  decay.  A cut-off  wave  as  the  one  with  ratio  0.7  is  also  well  simulated,  with  a linear  decay 
of  0.5237  dB/m  for  downstream  propagating  waves  and  0.5127  dB/m  for  upstream  propagating  waves 
( versus  a theoretical  decay  of  0.53 1 5 dB/m  ).  Nevertheless,  the  waves  with  % ratio  of  5.0  present  deficient 
prediction.  All  the  simulations  have  the  same  circumferential  wavelength  and  number  of  points  per 
circumferential  wavelength  ny  but  different  frequency  to  have  the  desired  cut-off  ratio.  This  change  in 
frequency  changes  the  propagation  characteristics  of  the  wave  and  therefore  a lower  number  of  point  per 
axial  wavelength  are  available  as  the  cut-off  ratio  increases  (for  inlet  perturbation). 


B ratio 

1.1 

3.0 

5.0 

Points  per  axial  wavelength  nx 

3569 

40 

22.5 

Table  1 Points  per  axial  wavelength  for  f,  =0. 7,  1.1,  3.0,  5.0  with  Nx=320 points 


Looking  at  Table  1 and  Figure  3 it  is  shown  that  the  accuracy  of  the  simulation  decreases  with  the  number 
of  points  per  axial  wavelength,  nx. 

In  addition,  this  accuracy  of  the  wave  propagation  prediction  shown  above  decreases  if  a coarser  mesh  is 
used  for  the  simulation,  with  different  worsening  depending  on  the  £,  ratio  (due  to  the  different 
propagation  angles  and  the  resulting  number  of  points  per  wavelength  ).  This  last  effect  is  clearly  seen  in 
Figure  4,  which  shows  the  results  of  upstream  and  downstream  running  waves  for  inlet  and  outlet 
perturbations  with  ^=3.0,  ny=20,  40  and  80  points  per  wavelength,  and  Nx/Ny=  1 .0 
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nrt  MOPiCiTION.  Cut  off  ratio  3.0.  HKSU  SeNSITlVITT 


Figure  4:  Mesh  sensitivity  in  the  simulation  of  the  propagation  of  a 3.0  cut  off  ratio  wave 


The  figures  shown  above  illustrate  how  the  numerical  results  of  the  wave  propagation  depend  strongly  on 
both  the  number  of  points  per  axial  and  circumferential  wavelength.  In  an  effort  to  unify  the  effects  of  ny 
and  nx,  Wilson  [ 9 ] proposed  the  analysis  of  wave  propagation  by  means  of  an  average  points  per 
wavelength 


f 


2,2 

x + ny 


(15) 


parameter  that  weights  the  least  defined  direction. 

The  analysis  of  all  the  simulations  computed  shows  that  the  dB  decay  per  wavelength  depends  strongly 
on  the  average  points  per  wavelength,  in  a correlation  close  to  the  3rd  order  between  the  error  and  the 
points  per  wavelength.  Taking  0.050  dB  decay  per  wavelength  as  threshold  for  acceptable  accuracy,  at 
least  40  average  points  are  required  for  the  accurate  simulation  of  propagation  waves.  From  this  average 
points  per  wavelength  needed,  it  raises  as  the  final  conclusion  that  at  least  50  points  per  wavelength  in 
circumferential  and  axial  directions  are  required  for  acceptable  accuracy. 

The  effect  of  the  mach  numbers  in  the  range  0.3<M<0.7  can  be  included  in  the  above  mentioned  analysis, 
as  the  damping  suffered  by  a wave  (either  upstream  or  downstream  potential  wave  as  well  as  vorticity 
waves)  traversing  the  physical  domain  is  moderated.  In  this  case,  the  effect  of  the  mach  number  will  be  a 
change  in  the  axial  propagation  characteristics  of  the  wave,  that  can  be  covered  with  the  mesh  analysis 
performed  above. 

However,  when  Mach  number  is  higher  than  0.7,  demonstrated  that  upstream  potential  fields  arc  highly 
damped  due  to  the  slow  propagation  velocity  of  such  waves.  In  such  a case,  the  previous  accuracy 
analysis  is  not  valid.  In  order  to  overcome  this  damping,  matricial  viscosity  was  implemented  in  the  code, 
to  reduce  damping  of  the  upstream  running  waves.  If  this  matricial  viscosity  is  used,  then  the  accuracy 
analysis  shown  before  is  also  valid  for  high  mach  numbers.  Results  of  these  simulations  with  the  MU  S T 
code  can  be  found  in  [ 6 ]. 


4 Flat  plate  response 

After  analysing  the  characteristic  parameters  of  the  wave  propagation  when  simulated  with  CFD  codes, 
the  next  step  is  the  interaction  of  a propagating  wave  with  an  isolated  profile. 

Before  going  into  detail  of  the  CFD  results,  an  analysis  of  the  parameters  involved  in  the  problem  is 
performed. 
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The  physical  characteristics  of  a flat  plate  response  to  unsteady  perturbations  depend  on  the  following 
variables 

® = f(lj,d,Ay,Ux,Uy,a>,P,T,Rg,r,U)  (16) 


that  involves  longitudes,  time,  mass  and  temperature.  Taking  the  pitch  distance  between  blade  rows  d,  the 
flow  velocity  U,  the  speed  of  sound  and  Rg  for  non  dimensionalization,  the  following  parameters  are 
found. 


® = f(-^,1-,^-,Mx,My,St,'Re)  (17) 
add 

fl 

where  St  stand  for  the  Strouhal  number  St  = - 

U 

Tyler  and  Sofrin  first  studied  the  interaction  between  blade  rows  [ 10  ],  showing  that  ratio  of  the 
wavelength  to  the  passage  pitch  can  only  take  certain  values  depending  on  the  blade  and  vane  number. 

The  theory  presented  in  section  2 can  be  extended  to  rotor-stator  interaction,  if  the  span  of  the  blade  rows 
is  small  compared  with  the  hub  radius,  which  is  an  standard  approximation  in  turbomachinery.  In  this 
case,  the  curvature  of  the  duct  passage  is  small,  and  therefore  circumferential  direction  can  be  assumed  to 
be  the  y direction,  with 


y = 

kyj 


R0 

- m j6 


( 18) 


Any  unsteady  variable  take  the  following  form 


The  wave  field  in  the  stator  frame  must  be  invariant  when  Bis  increased  the  angular  pitch  of  the  stator 
A=2tc  N and  t is  increased  by  A /£2,  due  to  the  unsteady  swirling  perturbations  with  angular  velocity  fi  . 

This  forces  the  mode  number,  m,  to  be 


m = kV  —nB  (20) 


In  the  stator  frame  of  reference,  the  response  should  be  circumferentialy  periodic  with  period  A=2jt  /V. 
Imposing  that  in  (19),  gives  the  following  expression  for  the  interblade  phase  angle  (IP A). 


V Xy 


Taking  this  into  account,  expression  (17)  takes  the  following  form 


® = f a,Mx, My, St,  Re)(22) 

d d 

As  done  in  section  3,  the  Reynolds  number  will  not  be  taken  into  account  as  the  study  will  be  limited  to 
potential  disturbances. 

Within  the  different  profile  characteristics,  the  flat  plates  present  the  advantage  that  a semi-analytical 
solution  of  the  unsteady  response  can  be  obtained  with  the  approach  presented  by  Whitehead  [ 5 ] by 
means  of  a linearized  solution  of  the  fluid  equations  (LINSUB).  Therefore,  this  study  will  make  use  of 

this  results,  and  the  effect  of  the  profile  — thickness  to  the  pitch  length  will  be  left  for  further  studies. 

d 
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The  effect  of  the  Strouhal  number  has  a quite  clear  physical  sense,  meaning  the  ratio  between  the 
characteristic  time  of  the  unsteadyness,  in  this  case  -y  of  the  incoming  wave,  and  the  characteristic  time 

for  the  fluid  to  travel  along  the  length  of  the  plate.  Results  from  LINSUB  showing  this  effect  on  the  plate 
response  are  shown  in  Figure  5. 


Ejfrct  oj  the  Strouhal  number.  M:  0.4.  d/l : 1.0 


Figure  5:  Effect  of  the  Strouhal  number  St  on  a flat  plate 


The  number  of  peaks  in  the  flat  plate  is  approximately  equal  to  the  Strouhal  number,  reflecting  that  the 
flow  changes  St  times  before  leaving  the  flat  plate  channel. 

However,  for  low  or  high  mach  numbers,  this  characteristic  is  distorted  due  to  the 
compressibility/incompressibility  effects.  Figure  6 illustrate  this  effect  for  perturbations  with  St=5.0  and 
IPA  270  degrees. 


Figure  6:  Effect  of  the  Mach  number  on  a flat  plate  response 


The  response  of  the  plate  at  high  mach  numbers  presents  more  peaks  but  the  peak  to  mean  level  of  the 
response  is  lower,  while  the  opposite  behaviour  takes  place  at  low  mach  numbers.  This  is  thought  to  be 
due  to  compressibility  effects  damping  the  response  while  making  easier  the  transmission  of  pressure 
waves  at  high  mach  numbers. 

In  ( 2 1 ) was  shown  that  the  IPA  depends  on  the  number  of  blades  and  vanes  of  the  stage.  Additionally 
the  IPA  also  reflects  the  ratio  of  the  circumferential  wavelength  of  the  incoming  disturbance  to  the 


(SYA)  4-9 


distance  between  the  flat  plates.  The  effect  of  this  ratio  is  shown  in  Figure  7,  for  the  case  of  a wave  with 
£=3.0,  and  a flat  plate  cascade  of  d/1  =1  and  a Strouhal  number  of  6.836 


Fiat  nwponst  fo  incoming  potential  diaturbanct 


Figure  7:  Flat  plate  response  for  different  interblade  phase  angles 


As  previously  explained,  with  M=0.4  the  amplitude  modulation  along  the  blade  row  presents  six  peaks 
for  all  the  interblade  phase  angle  analysed.  Nevertheless,  as  the  interblade  phase  angle  of  the  wave  is 
increased,  the  peak  to  mean  value  of  the  amplitude  modulation  increases,  to  a maximum  at  360  degrees.  If 
the  phase  profile  is  analysed,  it  can  be  seen  that  the  phase  modulation  is  similar  in  all  the  cases  studied  but 
in  the  case  of  360  degrees,  for  which  a different  trend  is  obtained. 

According  to  Woodley  & Peake  [ 1 1 ] [ 1 2 ],  when  unsteady  excitations  interact  with  a blade  row  (either 
vorticity  waves  convected  due  to  wakes  or  downstream-propagating  pressure  waves)  at  high  frequencies, 
resonance  on  the  blade  row  can  take  place  giving  a large  amplitude  response.  This  kind  of  behaviour  is 
known  as  Parker-mode-type  resonance,  and  according  to  [ 1 1 ] [ 12  ],  the  resonance  appears  only  under 
certain  combination  of  parameters  (Mean  flow  Mach  number,  stagger  angle,  frequency,  and  blade  and 
vane  number). 

To  test  the  capability  of  the  codes  to  simulate  the  response  of  the  flat  plate,  a circumferential  mode  30 
resulting  from  the  interaction  of  90  rotating  blades  with  60  staggered  flat  plates  was  computed.  The  mean 
flow  mach  number  of  0.4  giving  a £ ratio  of  3.0  with  and  interblade  phase  angle  of  180  degrees.  The 
Strouhal  number  is  set  to  6.836  and  the  pitch  to  chord  ratio  to  1 .0.  In  order  to  ensure  a good  propagation 
characteristics,  the  results  from  the  initial  investigation  are  used  to  set  the  grid  definition  upstream  and 
downstream  of  the  wave. 


Figure  8:  Pressure  field for  a 60 plates  configuration  with  I.P.A  180  ® M=0.4,  d/l=1.0  and  St=8.636 

In  Figure  8 the  resulting  flow  field  is  shown.  The  wave  propagates  with  an  angle  close  to  75°,  impinging 
in  the  flat  plate.  The  wave  keeps  on  propagating  in  the  channel  between  plates  and  is  radiated  downstream 
of  the  row  where  multiple  interaction  between  modes  gives  the  complex  pattern  of  the  figure. 
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Figure  9:  Flat  plate  response  for  incoming  disturbance  with  interblade  phase  angle  180  degrees.  Top: 
Actual  Values.  Bottom:  Percentage  Relative  Error 


The  results  of  the  flat  plate  response  are  shown  in  Figure  9. 

The  top  graphic  shows  how  the  interaction  between  the  flat  plates  and  the  wave  is  well  simulated  by  both 
the  linear  and  non  linear  version  of  the  MU2S2T  code.  The  six  peaks  resulting  on  the  response  are 
reproduced  by  all  the  computations,  and  good  agreement  with  the  phase  is  achieved.  The  greater 
differences  take  place  close  to  the  leading  edge,  probably  due  to  the  difficulties  on  the  numerical 
modelization  of  the  leading  edge.  Test  cases  with  lower  Strouhal  numbers  have  shown  lower  effects  of 
the  boundaries  on  the  plate  response.  In  this  case,  due  to  the  large  Strouhal  number,  the  flat  plate  response 
presents  high  complexity,  making  the  influence  of  the  leading  edge  significant.  In  the  rest  of  the  plate 
length,  differences  are  below  1 0 % of  the  values  predicted  by  LINSUB,  with  lower  values  as  the  mesh 
definition  is  refined.  In  Figure  9 is  also  reflected  the  higher  accuracy  of  the  linear  with  respect  to  the  non- 
linear version  of  the  code  for  the  same  number  of  points  in  the  flat  plate.  Table  2 displays  the  mean  error 
values  for  the  amplitude  and  phase  profiles  with  respect  to  the  LINSUB  value 


CASE 

Points  per 
Peak 

Amplitude 

Phase 

33 

6.3  % 

1 1 .5  % 

58 

4.6  % 

12.0% 

MUST  125  points 

20.8 

9.8  % 

20.9  % 

58 

7.3  % 

1 7.2  % 

Table  2:  Mean  relative  error  values  of  MUST  codes  with  respect  of  LINSUB  as  function  of  the  points  in 
the  flat  plate.  Interblade  Phase  Angle  180  deg 


A more  complex  test  case  arises  from  the  same  perturbation  with  an  interface  blade  angle  of  360  degrees, 
as  presented  in  Figure  7.  In  this  case,  the  resonance  makes  the  peak  to  mean  response  amplitude  larger, 
while  keeping  the  number  of  peaks  constant  with  respect  to  the  test  case  previously  shown. 
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An  interblade  phase  angle  of  360°  comes  from  a circumferential  mode  -30  resulting  from  the  interaction 
of  a 90  blade  configuration  with  90  flat  plates.  The  Mach  number  0.4,  with  Strouhal  number  of  6.836  and 
d/c=l  .0  are  kept  as  in  the  previous  test  case. 


Figure  10:  Pressure  field  for  a 90  plates  configuration  with  I.P.A  360,  M:  0.4,  d/l:  1.0  and  St:  8.636 

Figure  10  presents  the  resulting  flow  obtained.  As  in  the  previous  test  case,  the  potential  wave  impinges 
in  the  flat  plate  row  with  an  angle  close  to  75  degrees.  Due  to  the  higher  pitch  distance  between  flat 
plates,  the  reflection  pattern  in  the  channel  in  between  the  plates  is  simpler  than  in  the  180°  IPA,  with  the 
travelling  wave  pattern  still  visible.  Radiation  downstream  of  the  flat  plate  row  gives  a complex  field  of 
interactions. 


Figure  1 1 : Flat  plate  response  for  incoming  disturbance  with  interblade  phase  angle  360  degrees.  Top: 
Actual  Values.  Bottom:  Percentage  Relative  Error 
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The  Figure  1 1 shows  the  flat  plate  response  compared  with  the  results  from  LINSUB  for  the  linear  and 
the  non  linear  codes. 

The  top  graphic  shows  how  the  interaction  between  the  flat  plates  and  the  wave  is  well  simulated  by  both 

2 2 

the  linear  and  non  linear  version  of  the  MU  S T code,  accurately  predicting  the  six  peaks  of  the 
amplitude  response.  The  differences  reported  in  the  area  close  to  the  leading  edge  also  appear  in  all  these 
simulations.  Nevertheless,  the  correct  peak  to  mean  value  is  not  reached  with  none  of  the  simulations,  not 
even  with  a finer  mesh  of  500  points  in  the  flat  plate. 

The  phase  profile  presents  great  differences  in  the  area  of  the  trailing  edge,  where  none  of  the  simulations 
is  capable  to  reproduce  the  results  obtained  with  LINSUB.  The  change  in  the  trend  of  the  phase  variation 
is  a phenomena  observed  already  in  other  simulations,  but  that  effect  was  removed  with  mesh  refinement. 
Nevertheless,  in  this  test  case,  most  likely  due  to  the  high  complexity  of  the  plate  response,  the  mesh 
refinement  did  not  help  to  any  improvement  on  the  predictions  accuracy. 

The  relative  errors  with  respect  to  the  LINSUB  simulation  clearly  show  how  the  differences  along  the 
plate  are  low  (in  the  order  of  20-25  %),  except  in  the  local  maxima  and  minima  of  the  amplitude  response, 
where  the  relative  error  presents  much  larger  peaks. 

Table  3 displays  the  mean  error  values  for  the  amplitude  and  phase  profiles  with  respect  to  the  LINSUB 
value 


CASE 

Points  per 

Peak 

Amplitude 

Phase 

■KUBKESSa 

25 

26.2  % 

59.8% 

MUSTL  350  points 

58 

13.9% 

59.1  % 

MUSTL  500  points 

83 

15.5% 

60.1  % 

MUST  1 50  points 

25 

28.4  % 

66.3  % 

MUST  350  points 

58 

26.8  % 

65.4  % 

Table  3 : Mean  relative  error  values  of  MUST  codes  with  respect  of  LINSUB  as  function  of  the  points  in 
the  flat  plate.  Interblade  Phase  Angle  ISO  deg. 

It  can  be  seen  how  in  both  the  linear  and  non  linear  versions  of  the  MU2S2T  code,  the  relative  error 
decreases  when  the  mesh  is  refined  from  150  to  350  points.  Surprisingly,  when  a very  refined  mesh  with 
500  points  is  employed,  a worsening  in  the  solution  is  found. 

In  all  the  simulations  both  amplitude  and  phase  differences  respect  to  LINSUB  predictions  are  thought  to 
be  due  to  leading  and  trailing  edge  simulations.  The  plate  edges  effects  on  the  plate  response  have  been 
pointed  out  in  [ 1 1 ] [ 12  ] to  play  an  important  role  in  the  resonance  behavior  of  cascades  systems. 


5 Conclusions 


The  basic  principles  of  the  wave  propagation  and  flat  plate  interactions  simulation  have  been  shown. 

It  has  been  demonstrated  that  both  cut-on  and  cut-off  waves  can  be  simulated  with  CFD  codes.  In  order  to 
ensure  a correct  modelization  of  the  problem,  fully  unstructured  grids  with  2D  non  reflecting  boundary 
conditions  are  recommended,  since  they  have  shown  to  have  a major  impact  on  the  accuracy  of  the 
unsteady  flow  simulation.  The  main  parameter  of  the  accurate  wave  propagation  simulation  is  the  number 
of  points  per  wavelength,  both  in  circumferential  and  in  axial  direction.  It  has  been  shown  that  the  effect 
of  Mach  number,  flow  angle,  frequency  and  circumferential  mode  can  be  accurately  predicted  defining 
sufficient  grid  points  per  circumferential  and  axial  wavelength.  In  case  of  high  Mach  numbers,  this 
definition  could  be  not  enough  to  give  the  desired  accuracy,  and  matricial  viscosity  should  be  used. 

In  the  case  of  unsteady  response  of  flat  plates,  the  driving  parameters  have  been  shown.  For  the  range  of 
Mach  numbers  characteristic  of  low  pressure  turbines,  the  number  of  amplitude  peaks  can  be  determined 
with  the  Strouhal  number.  The  amplitude  of  such  response  have  be  shown  to  be  dependent  on  the 
interblade  phase  angle  of  the  wave.  Nevertheless,  60  points  per  response  peak  have  been  shown  to  give  a 
sufficiently  accurate  response  for  all  the  interblade  phase  angles. 
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The  behaviour  of  a non  linear  unstructured  Navier  Stokes  solver  for  simulating  wave  propagation  and 
transmission  through  a flat  plate  cascade  have  been  studied  and  the  results  have  been  shown  to  reproduce 
the  same  level  of  accuracy  than  the  linear  version  of  the  code  (but  with  some  more  mesh  refinement 
needed).  The  necessity  of  the  implementation  of  2D  UNRBC  for  minimize  the  reflections  of  the  wave  has 
also  been  demonstrated. 
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7 Nomenclature 

c:  Speed  of  Sound,  m/s 
d:  Pitch  distance,  m 
k:  Wave  number 

kn,j:  Axial  eigenvalue,  axial  wave  number,  m ' 
kyj:  Circumferential  eigenvalue,  circumferential  wave  number,  m ' 
m:  Circumferential  mode 
n:  Time  harmonic 

nx:  Number  of  points  per  axial  wavelength. 

ny:  Number  of  points  per  circumferential  wavelength. 

nav:  Average  points  per  wavelength. 

1:  Profile  length,  m 
t:  Profile  thickness,  m 
B:  Blade  number 

CFL:  Courant  Friedrichs  Lewy  parameter 
M:  Mach  number 

Nx:  Number  of  points  in  axial  direction 

Ny  : Number  of  points  in  circumferential  direction 

R:  Radius,  m 

Re:  Reynolds  Number 

Rg:  Gas  constant 

St:  Strouhal  number 

T:  Temperature,  K 

U:  Mean  flow  velocity,  m/s 

V:  Vane  number 

2 2 

MU  -S  T:  Multirow  Unsteady  Unstructured  Specific  Solver  for  Turbomachinery 

(D:  Angular  frequency,  rad/s 

(Or:  Resonance  frequency,  rad/s 

Xx:  Axial  wavelength,  m 

Xy:  Circumferential  wavelength,  m 

o:  Interface  Blade  Angle  (IP A),  deg 

£,:  Cut-off  ratio 

p:  Density,  kg/m3 

X-  Gas  Constant 

ji:  Viscosity  coefficient 

a2o4:  Artificial  viscosity  coefficients 

Q:  Shaft  rotational  speed,  rad/s 

4/:  Unsteady  variable 

A:  Angular  pitch,  deg. 


(SYA)  4-14 


References 


[ 1 ] A.  G.  Wilson.  Application  of  CFD  to  Wake/Aerofoil  Interaction  Noise-  A Flat  Plate  Validation  Case. 
AIAA  2001-2135.  May  2001. 


[ 2 ] F.  Kennepohl,  G.Kahl,  K.  Heinig.  Turbine  Blade/Vane  Interaction  Noise:  Calculation  with  a 3D 
Time-Linearised  Euler  Method.  AIAA-2001-2152.  May  2001. 


[ 3 ] B.  Lakshminarayana,  A.  Chernobrovkin,  D.  Ristic.  Unsteady  Viscous  Flow  Causing  Rotor-Stator 
Interaction  in  Turbines.  Journal  of  Propulsion  and  Power.  Vol.  16,  No,  September-October  2000. 


[ 4 ] G.A.Gerolymos,  D.  Vinteler,  R.  Haugeard,  G.  Tsanga,  I.  Vallet.  On  the  Computation  of  Unsteady 
Turbomachinery  Flows.  Part  2-  Rotor/Stator  Interaction  using  Euler  Equations.  AG  ARD  CP-57 1 . May 


1995. 


[ 5 ] D.S.  Whitehead.  Classical  Two  Dimensional  Methods.  Manual  on  Aeroelasticity  in  Axial  Flow 
Turbomachines , AGARD- AG-298,  Vol.  1,  Chap.  Ill,  1987 


[ 6 ] R.  Corral,  M.A.  Burgos,  A.  Garcia.  Influence  of  the  Artificial  dissipation  Model  on  the  Propagation 
of  Acoustic  and  Entropy  Waves.  45th  ASME  Gas  Turbine  and  Aeroengine  Congress.  8-1 1 May  2000 


[ 7 ] M.  B.  Giles.  Nonrefecting  Boundarv  Conditions  for  Euler  Equation  Calculations.  AIAA  Journal, 
Vol.  28.  No.  12,  pp  2050-2058 


[ 8 ] A.  Wilson.  A method  for  Deriving  Tone  Noise  Information  from  CFD  Calculations  on  the 
Aeroengine  Fan  Stage.  To  appear.  2000 


[ 9 ] A. Wilson.  TurboNoiseCFD  Deliverable  1.6:  Review  of  Generic  Studies.  15  April  2001 


[ 10  ] J.M.  Tyler,  T.G.  Sofrin.  Axial  Flow  Compressor  Noise  Studies.  SAE  Transactions,  Vol  70,  pp  309- 
332,  1962 


[11]  B.M. Woodley,  N.  Peake.  Resonant  acoustic  frequencies  of  a tandem  cascade.  Part  I.  Zero  relative 
motion.  J.  Fluid  Mech,  Vol  393,  pp  215-240,  1999 


[ 12  ] B.M.  Woodley,  N.  Peake.  Resonant  acoustic  frequencies  of  a tandem  cascade.  Part  2.  Rotating 
Blade  Rows.  J.  Fluid  Mech,  Vol  393,  pp  241-256,  1999 


(SYA)  4-15 


Reference  # of  Paper:  4 

Discusser’s  Name:  Dr.  Jan  Delfs 

Author’s  Name:  Mr.  Pedro  de  la  Calzada  Mazeres 

Question: 

What  is  the  reason  not  to  use  a CAA  approach  such  as  that  described  by 
Professor  Tam? 

Answer: 

We  are  primarily  interested  in  the  prediction  of  the  tonal  noise  due  to  stator-rotor 
interaction  in  turbines.  For  that,  we  need  to  solve  first  for  the  aerodynamics  in  order  to 
accurately  simulate  the  effect  of  both  the  potential  field  and  the  vortical  wake 
convection.  Our  intention  is  to  reduce  the  number  of  codes  used  in  the  simulation  of  the 
unsteady  aero  thermal  behavior  of  turbine  stages.  Therefore  we  are  concentrating  our 
efforts  to  have  only  one  code  with  the  necessary  features  implemented  to  predict  the 
unsteady  and  steady  flow  fields.  We  have  updated  our  CFD  code  for  the  steady  aero- 
thermal  field  in  turbines  and  we  want  to  keep  this  experience  and  improve  the  base  code 
further.  Therefore  we  are  working  on  the  code  in  the  direction  of  CAA  by  improving  the 
boundary  conditions  and  the  artificial  viscous  damping. 

Additionally,  the  modes  and  frequencies  in  which  we  are  interested  have  long 
wavelengths  on  the  order  of  the  pitch  in  the  circumferential  direction.  Thus,  the  mesh 
refinement  needed  to  resolve  the  viscous  boundary  layer  and  wakes  is  probably  finer 
than  that  needed  to  resolve  the  acoustics.  So,  we  are  happy  to  solve  for  the  unsteady 
behavior  with  the  same  mesh  and  code  we  use  for  the  steady  aero-thermal  field. 


Discusser’s  Name:  Mr.  Jose  M.  Riola  Rodriguez 

Author’s  Name:  Mr.  Pedro  de  la  Calzada  Mazeres 

Question: 

How  long  does  the  nonlinear  code  take  to  run  compared  to  the  linear  version? 
Answer: 

We  have  found  a factor  of  two  between  the  two  codes.  In  order  to  have  an  idea 
of  the  total  time  needed  for  the  flat  plate  cases  presented  here:  for  a mesh  with  around 
20,000  nodes,  the  linear  version  runs  overnight  and  the  nonlinear  version  needs  a full 
day. 


